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Abstract. We study a simple one-dimensional quantum system on a circle 
with n scale free point interactions. The spectrum of this system is discrete and 
expressible as a solution of an explicit secular equation. However, its statistical 
properties are nontrivial. The level spacing distribution between its neighbor- 
ing odd and even levels displays a surprising agreement with the prediction 
obtained for the Gaussian Orthogonal Ensemble of random matrices. 

1. Introduction 

One possible approach to study the chaotic properties of bounded time-independent 
quantum system is based on results of the Random matrix theory (RMT). Bohi- 
gas, Giannoni and Schmit pQ conjectured that the local fluctuations of quantum 
energy levels of such systems display universal properties. Following this conjecture 
the level spacing distribution of integrable systems with more than one degree of 
freedom is expected to be Poissonian, while the distribution for systems that are 
classically chaotic is described by the random matrix theory (see [5]). (Integrabil- 
ity in the quantum case means that the eigenvalues and eigenvectors of the system 
can be evaluated by solving a simple algebraic equation.) The specific distribution 
depends in this case on the symmetry of the system and is described by the Gauss- 
ian orthogonal, unitary or symplectic ensemble respectively. Nowadays there is an 
impressive amount of evidence for the relevance of this conjecture and so for the 
link between the random matrix theory and the quantum behaviour of classically 
chaotic systems. 

The assumption about the number of degrees of freedom is, however, quite im- 
portant since the conjecture doesn't hold in one dimension. All classical time 
independent one dimensional systems are integrable. Nevertheless there are quasi 
one-dimensional quantum systems (the quantum graphs) that display random ma- 
trix properties. Recent papers give a wide theoretical and experimental evidence for 
such behaviour in various interconnected quantum graphs [3 [H [6] . The quantum 
dynamic is here related to the nontrivial topological structure of the graph and the 
random matrix description of the spectral statistics of the quantum graph is based 
on this fact. 
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Quantum graphs are topologically nontrivial since each node is connected with 
at least three neighbors. This is why these systems are regarded as quasi one 
dimensional. In a very one dimensional system each point is topologically related 
only to its two neighbors: to the one on its right and to the one on its left. Our 
aim here is to construct a virtually one dimensional system (i.e. each node has 
only 2 neighbors) displaying a random matrix behavior. The work is based on the 
paper [llj . However there is a great difference in results. We focus on the weak 
coupling limit and show that in this case, the level-spacing spectral statistic display 
two fundamentally different types of behavior. The statistical properties of the odd 
level spacings are described by the random matrix theory while the even spacings 
do not show this behavior. 

We deal with a one dimensional particle on a circle with scale free point inter- 
actions localized at the points < X\ < . . . < x n < 2ir. The scale free point 
interaction ([7J [B]) localized at the origin is described by a one parameter family 
of boundary conditions: -/(0_) = /(0+), a/'(0_) = /'(0+). It has the amazing 
property that the corresponding reflection / transmission coefficients do not depend 
on the particle energy. The scale free point interaction is also expressible as a limit 
of short-range potentials, see [9l HOj. 

The Hamiltonian of the system is represented by the one dimensional Laplace 
operator 

with a set of boundary conditions: 

(2) /(2tt_) = /(0+), /'(2tt_) = /'(0+), 

(3) -f(xk-) = f{xk+), a/'(x k _) = f'(x k+ ), k = l,...,n, 
a 

where a 6 (0, +oo) is a parameter. We denote the domain of T a as D(T a ). The 
periodic boundary conditions ([2]) define a circle. The second pair ([3]) describes the 
point interactions. It can be easily shown (see the theory of self adjoint extensions 
e.g.[H]) that the operator T a is self adjoint for every a. Moreover, its spectrum 
is purely discrete containing eigenvalues of finite multiplicity without limit points. 
Note that for a = 1 the Hamiltonian describes simply a free particle on a circle. 
Furthermore in the limits a — > 0, a — > +oo the system decouples into n independent 
subsystems with the Dirichlet boundary condition on one side and the Neumann 
condition on the other side of the subintervals (x k , Xk+i)- 

2. The spectrum 

In what follows we will focus on the positive part of the spectrum denoting 
x a = and x n+ \ = 27r. An eigenvector ifk,k > corresponding to the energy 
E = k 2 writes on the subinterval (xj, Xj+i), j = 0, . . . , n as 

(4) <p k (x) - A 3 {ky kx + B,(k)e- tkx , 

with some coefficients Aj and Bj. The vector (fk(x) belongs to the domain D(T a ) 
if it satisfies the boundary conditions and ([3]) . This leads to a relation between 
the neighboring coefficients Aj— i, flj— i and Aj, Bj 
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for each j = 1, . . . , n, with Cj being a 2x2 matrix 
(5) C 



If (l + a) (L-a)e- i2kx 



j 



2\ (i - a) e ak *i (i + a) 



The dependence between A n , B n and An, Bq is given by the periodic boundary 
conditions ([2]) and leads to a relation 

(6) f^=C^ 



Bn \B, 



with matrix 

(7) C P = 



e 



-i2k7r 



Combining all these relations together we obtain that the coefficients A Q , B de- 
termine an eigenvector of T a iff 

(8) (C P C n • • • d - I) = 0. 

So the eigenvalues of T a are obtained as solutions of the equation 

(9) (p(k) = det(C P C n ---C 1 -I) = 0. 

Since the matrices in the above determinant are of size 2x2, we can easily 
determine the explicit form of ^ from © and ([7]) by induction. The secular 
equation becomes 

n 

(10) ( P {k) = cos(2fc7r) - (1 - /3 2 )^ + (3 2 ^ cos(2fc(x fc - xi + vr)) + /? 4 . . . = 

k<l 

where 

I -a 2 



(11) 



1 + a 2 



Note that (TTu| is a finite series and it ends when the number of indices in the last 
sum reaches the highest even number m with m < n. So the last term is 

n 

(12) /3 m ^2 cos(2k(x il - x i2 +... + x im _ 1 - x im + n)). 

ii<...<i m 

Since the equation (fiTJj) contains only even powers of (3 the positive spectrum of 
T a is invariant under the replacement a <-> I /a (it leads to the change f3 <-> —0). 
In the free case (a = 1;(3 = 0) the equation (TIT)]) simplifies to 

(13) CpW = cos(2fc7r) - 1 = 
and the solutions are just 

(14) h = i, l e N, 

with the eigenvalues being doubly degenerate (a consequence of the rotational sym- 
metry of the system). 
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3. Statistical properties of the spectrum 

We start with the free case omitting its non degenerate ground state. Further 
we suppose that all the points < X\ < . . . < x n < i n+ i = 27r are rationally 
independent, i.e. we suppose that the equation 

n+l 

(15) mm = 0, mi E 1 

j=i 

has a solution only for m, = for all i. 

The states are two times degenerate. Since we want to compare the fluctuation 
properties of the energy levels with the random matrix theory we will not work 
directly with the eigenvalues Ei = kf but with the values e; = 2|fc;|. The point 
is that the values e; have not to be unfolded and are of density 1. From now on 
whenever we deal with the level spacing statistics we have the "spectrum" e; in 
mind. 

The free case is trivial: e2i-i — en = 21,1 = 1,2, .... Defining the level spacing 
si = ei + i — ei we immediately see that the probability density P(s) consists simply 
of two delta peaks P{s) = (S(s) + 5(s — 2))/2. The peak at comes from the 
degeneracy of the eigenvalues. 

The question is what happens in the case a ^ 1, i.e. when the point interactions 
are switched on. The degeneracy vanishes and the two delta peaks turn into smooth 
distributions. 

To see what happens in that case let us assume a small change of the parameter 

a : 

a = 1 + 6 

with \S\ < 1. The relation (fTl]) leads to 

So the secular equation (|10[) can be solved perturbatively supposing the roots kj 
in the form 

(16) kj = j + pXj, j 6NU{0} 

with |/3 Xj | <C 1. Substituting it into the equation ( (fTu| ) and using the properties 
of goniometric functions we get 

1 — 1 + . . 

2 2 




+(3 2 ^cos{2j{x k -xi)) \ 1- 

k<l 
n 

-/3 2 ^2sin(2j(x k - xi)) (2f3Xj{x k - x\ + 7r) - . . .] 

k<l 

The absolute term drops out and we get (in the order of j3 2 ) 

n 

-2n 2 X 2 + | + cos(2j(x fe - xi)) = 0. 

k<l 
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The solutions Xj are given by 
1 



(17) Xf=± 




^cos(2jx fc ) + ^ sin(2jx fc 



The degenerate eigenvalues split and become 

(18) k 2j ^=j + \l3\XJ, k 2j =j + \P\Xj 3 = 1,2,... 
So the unfolded level spacings are 

(19) S2J _ 1 =4|/3|A+ S2J =2-2|/3|(A+ +1 +A+). 

The odd spacings S2j-x contribute to the widening of the S~ peak localized at 
while the even spacings s 2 j broaden the 8 peak localized at 2. We will now show 
that the peak at is not only broaden but it acquires a shape prescribed by the 
random matrix theory. 

The positions x^ and n are altogether rationally independent. Hence Xi/n are 
also rationally independent irrational numbers. The periodicity of cosine leads to 
the relation 

cos(2jxi) = cos(2jxi mod 2tt) = cos(27r((jxi/7r) mod 1)), 

where mod means the remainder after integer division. From the theory of dis- 
tribution of sequences modulo one (see |13j ) we know that when j changes the 
arguments of cosines and sines in (|17|) behave like independent random variables 
uniformly distributed in [0, 27r). From this fact follows that the distribution of the 
sequence {cos(2 jxi)}^^ has a probability density function 

xe (-1,1), 



(20) /(*)= ^ 



0, I6l\(-l,l). 

Moreover for different positions Xi, i — 1, ...,n these sequences behave like in- 
dependent random variables. The same holds for {sm(2jxk)}j. So the sums 
*Efe=i cos (2jxk)}j, {Sfe=i sm (2j^fc)}j normalized by the factor V2/ ' \fn converge 
with increasing n to two normally distributed random variables. 

In fact there two variables are statistically independent. To see this we sketch 
the proof of this fact. It is based on the observation that the pair (cos(X), sin(X)), 
with X being uniformly distributed along (0, 2ir), can be equivalently represented 
in the form 



sgn(cos X) | cos X \ , sgn(sin X) \J\ — \ cosX| 

As one can easily prove, random variables sgn(cosX), sgn(sinX) and |cosX| are 
stochastically independent. So the statistical dependence of the cosine and sine 
is fully contained in their absolute values, while their signs are statistically inde- 
pendent. This observation together with the characteristic function technique used 
in the standard proof of the Central limit theorem shows finally the statistical 
independence of the above sums. 

Summarizing we see that in the weak perturbation limit the distribution of the 
odd spacings S2j-i is given by a square root of a sum of squares of two independent 
and normally distributed random variables. Properly normalized it is nothing but 
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the Wigner distribution P\y approximating the spacing distribution Pqo e of the 
Gaussian orthogonal ensemble of random matrices (GOE): 

(21) P w (s) = | se -^ 2 / 4 . 

Remember that during the limiting procedure the perturbation condition \/3Xj | <C 
1 has to be fulfilled for every n. Together with the structure of the formula (fl~7|) it 
leads to a condition |/3| <C 1/n. Since the exact Wigner distribution is obtained in 
the limit n — > oo the parameter [3 has to decrease to 0. It is an open question how 
this distribution behaves for increasing n and fixed (3. We suppose that for fixed j3 
the difference to the Wigner distribution first decreases with increasing n. Than, 
for some value of n, it reaches a minimum and starts to increase when n violates 
the condition <§; 1/n. 

The even spacings S2j behave differently. Substracting the constant 2 the relation 
(|19[) contains a sum of two variables with a Wigner distribution. So the spacing 
probability density behaves as ~ s 3 for small s and hence in a way that is not related 
to the random matrix theory. The next section shows that the distribution of the 
odd spacings follows the random matrix prediction even in the non perturbative 
regime. 

The topology of the system (a circle) is of fundamental importance. To show 
this we investigate an analogous system on a line segment. All the parameters of 
the system remain unchanged. The only difference is that the periodic boundary 
conditions ((2]) are replaced by the Dirichlet ones: /(0) = /(27r) = 0. The spectrum 
is again purely discrete and the eigenvalues are determined by the roots of the 
secular equation 

n 

Q L {k) = sin(2fc7r) + (3 sin(2k(x k - tt)) + 

fc 

n n 

(22) +(3 2 ^2cos{2k(x k -xi +tt)) + /3 3 ••• = 

k<l k<Km 

The free case (/3 = 0) leads to unfolded energies e/ = 1,1 = 1, 2, ... and the proba- 
bility density P(s) of the level spacings s/ = ei+\ — e; is just P(s) = S(s — 1). 

For weak coupling one applies the perturbation methods. With the assumption 
kj = j/2 + /3"fj + 0(f3 2 ), j e N the result is given by 

1 - 

(23) lj = 2Zz2sin{jx k )- 

fc 

The distribution of jj tends (after rescaling by the factor to a normal 

distribution with zero mean as n — > oo. So the unfolded level spacings, Sj — 
1 + 2/3(7^+1 — 7j), are in the weak coupling limit normally distributed and not 
related to the random matrix theory. 



4. Results 

We solved the equation © finding its first 10 6 roots for different values of a and 
n. The interaction positions have been chosen as 

(24) Xk = ^±± t k = l,2,..n 
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with pk denoting the k— th prime number. This choice is obviously rationally in- 
dependent and (fl"5|) holds. The presented results are, however, not influenced by 
this particular choice. For different positions satisfying (|15p the results remain 
unchanged. 

In the perturbative regime the odd spacings follow the Wigner distribution 
whereas the even spacings not. Interestingly enough the spacing statistics behave 
similarly also beyond the weak coupling limit. 

To compare the obtained level spacing distributions with the results known from 
the random matrix theory we have evaluated the difference functions 

+00 

(25) AF W = J (F(s)- F w (s)) 2 ds 




-t-00 



(26) AFgoe = J (F( s )-F G oE{s)) 2 ds 



where Fyy and Fqoe is the integrated spacing distributions for the Wigner and 
the exact GOE result respectively: 

s s 

(27) F w (s) = J P w {u)du- Fgoe(s) = J P G OE{u)du. 


As an approximation of the exact distribution Pgoe we used the Taylor ex- 
pansion up to the power 42 for small spacings and the Dyson asymptotic result 
otherwise. This approach is described in [2] chapter 4.9. It leads to a difference 
between the Wigner surmise and the GOE result: 

+00 

(28) A(F GOE -F w )= J (F GO E(s)-F w ( s )) 2 d s = 3.9280 • HT 5 . 



The distribution F(s) of the odd level spacings {52^-1}^ is defined in a usual 
way, 

1 N 

(29) F(«) = 

3=1 

with O being the Heaviside step function and N the number of spacings taken into 
account. Its derivative (in the limit TV — > 00) is the level spacing density. 

A typical result is shown on the figure [T] The solid line marks the Wigner 
probability density, (f2"Tj) . Stars show the a the level-spacing probability density for 
GOE as published in [14], Table A. 15. 

The dependence of the differences AFw and AFgoe on the parameter a is 
shown in the figure [2] 

It is not a surprise that for the weak coupling a ~ 1 the sequence Saj-i follows 
the Wigner distribution (|2"l"j) more closely than Mehta-Gaudin (GOE) distribution. 
This is just the consequence of the arguments presented in the previous section. It 
remains true, however, also for a up to ~ 1.35 i.e. far outside the weak coupling 
limit. Even more: there is a distinct minimum of AFw for a w 1.275. With 
increasing a the distribution starts, however, approach closer the GOE result then 
its Wigner approximation. For a w 1.4 the difference AFgoe displays a clear 
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Figure 1. The odd level spacing distribution of S2.7-1 for 47 point 
interactions with a = 1.001 and N = 1 • 10 5 is compared with the 
Wigner distribution (solid line) and the exact GOE result (stars). 




1.001 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 



Figure 2. Plotted are the differences AF W and AF GO e for 24 
point interactions and N — 10 6 . The dashed horizontal line repre- 
sents the difference between true GOE and Wigner - (|28|) . 

minimum with Agoe = 3.6457 • 10~ 7 . For this coupling the spacing distribution 
of the model very close to the true GOE statistic. For a > 1.5 the agreement 
however deteriorates step by step. This is understandable since the system with 
many non-periodically placed interactions is influenced by localization effects for 
strong enough coupling. For a — > oo the spectrum becomes a sum of n mutually 
embedded equidistant sequences and becomes Poissonian for rationally independent 
points Xk - see [T3] . 

On the other hand the distribution of even spacings S2j never becomes random 
matrix like, see figure [3l 

The level spacing statistic depends only on the neighbouring levels. More sensi- 
tive statistical properties, like for instance the number variance, describe the cor- 
relation of larger level sets. In our case, such large correlations are dominated by 
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Figure 3. The even level spacing probability density for 9 point 
interactions with a = 1.9 and N = 1 • 10 5 is compared with the 
Wigner and the GOE distributions. 

the simple geometric structure of the model and by the related periodic orbits (see 
|16j). This leads finally to oscillatory pattern like the one observed in the number 
variance and plotted on the figure [4] 




s 



Figure 4. The number variance [12 = (n — (n)) 2 evaluated for 
24 point interactions with a = 1.4 and N = 1 • 10 6 . The solid 
and dashed curves correspond to the predictions for the Gaussian 
orthogonal (GOE) and the Gaussian unitary (GUE) ensembles re- 
spectively. 

The spectral properties of the model are closely related to its topology. Remem- 
ber that the odd spacings originate from the splitting of the two times degenerate 
eigenvalues for the free case and are in such a way linked to the fact that the sys- 
tem has a topology of a circle. If we change the topology and replace the periodic 
boundary condition ([2|) with the Dirichlet ones /(0) = /(27r) = (i.e. we work on a 
line segment instead of on a circle) the spacing distribution changes fundamentally 
and never follows the random matrix theory. In this case the distribution of all 
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spacings (odd and even) is the same. Their common distribution is displayed on 
the figure [5] 

1 1 1 1 1 1 1 1 1 1 



0.8 - 




Figure 5 . The level spacing distribution of Sj for 9 point interac- 
tions on a line segment with a — 1.8 and N = 1 • 10 5 is compared 
with the Wigner and GOE distributions. 

To summarize: we have discussed a simple one dimensional quantum model 
that displays a level spacing statistics being in a surprising agreement with the 
predictions of the Gaussian orthogonal ensemble of random matrices. We show 
that this agreement is related to the topology of the system (circle) and vanishes 
when the topology is changed (circle is replaced by a line segment). 
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